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Abstract 

A high precision method to determine fnsion in the captnre of a-particles 
by nnclei is presented. For a-captnre by "^^Ca and snch an approach 

gives (1) the parameters of the a-nnclens potential and (2) fnsion probabil¬ 
ities. This method fonnd new parametrization and fusion probabilities and 
decreased the error by 41.72 times for a Ca and 34.06 times for a Ca 
in a description of experimental data in comparison with existing results. We 
show that the sharp angular momentum cutoff proposed by Glas and Mosel 
is a rough approximation, Wong’s formula and the Hill-Wheeler approach 
determine the penetrability of the barrier without a correct consideration of 
the barrier shape, and the WKB approach gives reduced fusion probabilities. 
Based on our fusion probability formula, we explain the difference between 
experimental cross-sections for a+^^Ca. and a-l-"^^Ca, which is connected with 
the theory of coexistence of the spherical and deformed shapes in the ground 
state for nuclei near the neutron magic shell N = 20. To provide deeper 
insight into the physics of nuclei with the new magic number N = 26, the 
cross-section for a Ca is predicted for future experimental tests. The role 
of nuclear deformations in calculations of the fusion probabilities is analyzed. 
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1. Introduction 


Understanding fusion in nuclear reactions is a fundamental problem of 
physics [H. A topic where such processes are important is the synthesis of 
superheavy nuclei Bi- Here, fusion has a crucial role, but its analysis is 
complicated. Deep insight into the physics of fusion can be obtained from 
the capture of the o-particles by the nuclei. 

Information about fusion in a-capture forms our understanding of inter¬ 
actions between the a-particles and nuclei at distances, where an important 
role is attributed to the formation of the nucleus from two colliding nu¬ 
clear objects inside the spatial nuclear region, a-nucleus interactions have 
been studied most intensively in the context of the a-decay of nuclei (see 


iJ3, 


various microscopic 


models 


macro- 


31-44 1, hssion models B 46|), and the scattering of 


experimental information 
scopic cluster models 
the a-particles off nuclei 47H51|. The physics of the fusion processes dur¬ 
ing a-capture have been investigated less deeply laMi. Evaluations of the 
a-particle capture rates show that they play an important role in nuclear 


reactions in stars 50, 52 


A prevailing approach for the determination of the capture cross-sections 
of the a-particle by the nucleus is based on calculations of the penetra¬ 
bility of the potential barriers, where information about the fusion is not 
included. While many approaches have been developed for the calculation 
of the penetrability, there is no generally accepted method to estimate the 
fusion. Experimentally, these reactions have not been studied deeply: we 
have the cross-sections for capture of the a-particles by nuclei ^°Ca, ^"^Ca 


531, ^^Co 


208 


Pb [^, and ^°®Be 55 


We shall be interesting in information about fusion that can be extracted 
from the available experimental data. Thus, the aim of this paper is to con¬ 
struct a new approach for obtaining such information. The idea of such an 


approach can be similar to the inverse theory of scattering jSg . Eberhard et 


al proposed a relation that gives information about fusion in the a-capture 
and compares calculated cross-sections with experimental data at selected 
energies |53|. Analyzing the existed experimental data, we observe an inter¬ 
esting difference in fusion between a-captures by the ^°Ca and “^^Ca nuclei, 
while the inclusion of the other nuclei into the analysis does not result in 
principally new questions. 
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A crucial point in such a task is the determination of the barrier penetra¬ 
bility. It turns out that even small improvements in the Wentzel-Kramers- 
Brillouin (WKB) formula of the penetrability require the reconsideration of 
the connection between boundary conditions and even the inclusion of the 
initial (or hnal) condition into the stationary picture of the studied process. 
After such a modihcation, the penetrability becomes more sensitive to the 
shape of the barrier, and it is already dependent on the potential outside the 
tunneling region. A fully quantum realization of the approach leads to the 
essential influence of the calculated penetrability on additional parameters 
for details, proofs, demonstrations). The role of such additional 


see 


parameters in the determination of the penetrability is larger in some orders 
in comparison with changes of the resulting probability as a result of any 
variations in the nuclear deformations and parameters of the interacting po¬ 
tentials if the penetrability is not calculated in the fully quantum approach. 
In this paper, we shall therefore restrict ourselves to the spherically symmet¬ 
ric a-nucleus potentials. However, the role of the nuclear deformations in the 
determination of the fusion probabilities in the a-capture task is analyzed 
Appendix F 


m 


We generalize the method of multiple internal reflections (MIR method. 


see |57H62|I to describe the capture of the a-particles by nuclei. We analyze 
the sharp angular momentum cutoff approach previously introduced by Glas 
and Mosel in j^, which has been widely used (for example, see 69]). 


We analyze the WKB approach, Wong’s formula [6^ and the Hill-Wheeler 
approach in calculations of the cross-sections. 

We focus on understanding reasons for difference between the cross-sections 
for a Ca and a Ca. The inclusion of the fusion process into the a- 
capture task allows us to essentially improve the agreement between the 
calculations and experimental data, and it provides a deeper picture of the 
a-capture than the sharp angular momentum cut off approach. We show 
that the MIR approach, in contrast to other approaches, satishes all tests of 
quantum mechanics, resolves the problem of indehnite additional parameters 
described above (hxing the result and excluding errors), and it is the most 
effective and accurate method for studying fusion in this task. We hnd a 
new formula for fusion probability and propose a-capture estimations for the 
^®Ca nucleus with the new neutron magic number N = 26 


67 
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2. Method 


2.1. Cross-section of a-capture and sharp angular momentum cut-off 

The cross-section of capture that includes the fusion of the a-particles 
with nuclei can be dehned as 53 


+°° 

acapture(i?) = —— E(2^ + 1) ^1 Pi, (1) 

2m E 

where E is the kinetic energy of the a-particle incident on the nucleus in 
the laboratory frame, E is the kinetic energy of the relative motion of the 
a-particle of the nucleus in the center-of-mass frame (we use T' ~ i?), m 
is the reduced mass of the a-particle and nucleus, Pi is the probability for 
fusion of the a-particle and nucleus, and Ti is the penetrability of the barrier. 
In the WKB approach, this coefficient is dehned as 

Twkb{E) = exp |-2y - ld(r)) drj, (2) 

R 2 


where R 2 and are the second and third turning points determining the 


10 - 


tunneling region! 

We shall consider the capture when the fragment has passed to the in¬ 
ternal potential well after tunneling through the barrier. Solving a set of 
classical motion equations, one can hnd the critical angular momentum 
when all trajectories at I < Cr lead to capture. Glas and Mosel introduced 


the sharp angular momentum cutoff [63|, l6^ , where 


Pi = 


1 at 
0 at I > I 


I C: lev, 


( 3 ) 


In some papers (for example, see 53|), an assumption that Ti = 1 for I < /gr 
is imposed (where /gr > /„), which transforms formula (|2]) into 


Cr capture (-/^ ) 


7rr 

2mE 


(2/cr -l- 1)^ 


( 4 ) 


'^The first turning point is defined as the intersection of the potential and energy near 
zero, where the centrifugal component is the main contribution 
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The critical momenta /cr are found after comparison of the cross-sections 


obtained by this formula with experimental data. Eberhard et al in [53 


obtained = 9.9 for and = 10.9 for a-f-'^^Ca. We shall reconsider 

this process for the quantum consideration that leads to a set of tunneling 
transitions at each allowed orbital momentum I according to formula ([T]). 

2.2. Multiple internal reflections method in the barrier penetrability determi¬ 
nation of an arbitrary shape 

We shall study the capture of the a-particle by the nucleus in a spherically 
symmetric scenario (the role of nuclear deformations in the determination of 


the fusion probabilities in the a-capture task is analyzed in Appendix F). 
To apply the idea of multiple internal reflections to study packet tunneling 
through complicated realistic barriers, let us consider the radial barrier of an 
arbitrary shape, which has successfully been approximated by a sufficiently 
large number N of rectangular steps: 


V{r) = 


Vl at Tmin < r < ri 
V 2 at ri < r < r 2 

Vjv at rjY_i < r < r„ 


(region 1), 
(region 2), 

(region N), 


( 5 ) 


where Vj are constants (j = l... A^jl. Let us denote the hrst region with a 
left boundary at point rmin (we denote it also as Tcapture), and we shall assume 
that the capture of the a-particle by the nucleus in this region occurs after 
its tunneling through the barrier. We shall be interested in solutions for the 
above barrier energies, while the solution for tunneling could be subsequently 
obtained by changing —)■ ki. A general solution of the wave function (up 
to its normalization) has the following form: 




( 6 ) 


^We used this approximation for barriers of the proton and a-decay for the number 
of nuclei, where the width of each step is equal to 0.01 fm, and we demonstrated the 


stability of the calculations of all amplitudes of the wave function and penetrability 57l |. 
We thus have effective tools for a detailed study of the quantum processes of tunneling 
and penetrability. 
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^ikir 


^—ikir 


X{r) = < 


ai + /?! e 

a^e^k 2 r ^ ^ 


at r^in < r < ri 
at ri < r < r 2 


(region 1), 
(region 2), 


CTtv—1 6 

p-ikffr 




+ e 


+ f^N-i e 

ikj^r 


—ikj^—ir 


at rN -2 <r < tat-i 
at Ttv-i < r < r^ax 


(region — 1), 
(region N), 

m 

where aj and [3j are nnknown amplitndes, Aj- and are nnknown am¬ 
plitudes of transition and reflection, Yi^iO, cp) is the spherical function, and 


kj = ^y2m{E — Vj) are complex wave numbers. We have hxed the normal¬ 
ization so that the modulus of amplitude of the starting wave equals 

unity. We shall search a solution to this problem by the multiple internal 
reflections approach. 

According to the multiple internal reflections method, the scattering of 
a particle on the barrier is sequentially considered by steps of propagation 
of the wave packet relative to each boundary of the barrier (the idea of this 
approach can be understood most clearly in the problem of tunneling through 


the simplest rectangular barrier, see |59l. l6ll. l62| . where one can hnd proof 


of this fully quantum exactly solvable method and analyze its properties). 
Each step in such a consideration of packet propagation is similar to one of 
the hrst independent 2N — 1 steps. From the analysis of these steps, we hnd 
recurrent relations for the calculation of the unknown amplitudes A^^\ ^r\ 
and for an arbitrary step with number n (here, index j corresponds 
to the number of region V)-, and the logic of the dehnition of these amplitudes 
can be found in Appendix A in j^). In the summation of these relations 
from each step, we impose the continuity condition for the full (summarized) 
wave function and its derivative relative to the corresponding boundary. 

According to the analysis of waves propagating in a region with an ar¬ 
bitrary number j on the arbitrary step, each wave can be represented as a 
multiplication of the exponential factor ^ and constant amplitude. In 
practical calculations, the difficulty consists in the determination of these 
unknown amplitudes. However, to make such calculations as easy as pos¬ 
sible for an arbitrary step, one can rewrite the amplitude of the wave that 
has transmitted through the boundary with number j as the product of the 
amplitude of the corresponding wave incident on this boundary and the new 
factor T'^ (i. e., the amplitude of transition through the boundary with num¬ 
ber j). The bottom index indicates the number of the boundary, and the 
upper sign or ” is the direction of the incident wave to the right 
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or left, respectively. We associate the amplitude of the reflected wave from 
the boundary with number j with the amplitude of the wave incident on 


this boundary via new factors Rf. The coefficients Tj^, T^.. .and i?f, 
Rf, Rf... can be found from the recurrence relations shown above We 




n± 


3± 


Tt 

rp± 

I -1-2 • 

rp± 

■ ■ ■ -‘-N- 

^ and Rf, 

Rf 

...R 

Tf 


2k, 


"j 

tt = 

3 

kj 

+ fcj+i 



3 

Rf 

kj 

- kj+i 

2ikjrj 


Rj = 

3 

kj 

+ kj^i 



3 


N-l 


as 




kj + fcj+i 
fcj+i — kj 
kj + fcj+i 


( 8 ) 


2_ ^-2ikj+irj 


Now, we consider the wave propagating in region with number j, which 
is incident from the right on the potential barrier with the right boundary 
at point rj_i (and left boundary at point ri). Let us hnd the wave reflected 
from this complicated barrier. This wave should combine all waves formed as 
a result of multiple internal reflections and propagations relative boundaries 
ri ... rj_i and leave such a barrier. We dehne the reflection amplitude Rf^i 
of such a summarized wave as 


= Rl_, + r,+.fi+T,-i(l + E (RtRj-ir) = RU + 


RJR-_, 


(9) 


Correspondingly, we also dehne 


R 


i+i 


= R 


~ / “1“00 ~ \ 

Tj_^_lRj {Rf+iRj = Rj+1 

'■ 'm — 1 ' 


-j+1 


T+ — 

^j+i - 


i+i i i+i' 

H-oo 


+ 


m=l 


1 - Rf+iR~' 


T.{RURjr) 

^ m=l ' 


T. 




, +00 

TrT-_Ai+ E 

^ m=l 


'3 


1 - RAiRj 

1 - R-_^Rf' 


( 10 ) 

In such a summation, we have recurrent relations, which connect all ampli¬ 
tudes. We now choose the following values 

Rfr-l = -RjV-l) -^1 = -^1 ) ) ^AT-l = ^AT-l) (H) 

as our starting point and consequently calculate all amplitudes Rfi _2 ■ ■ ■ Rf, 
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R2 ■ ■ ■ Rn-1 T2'" ... We find the coefficients aj and f 3 j: 


+ 00 


■|“LXJ / \ ~ / ~1“00 ~ \ 

E «r = ^- 1(1 + s (RtRj-i)"") 

1=1 ^ m=l ' 

H“00 / \ yw / ~1“00 yw \ 

= E 1(1 + E 

n=l ^ m=l ^ 


aj — 


^. 7-1 




T+ 

] 

rp+ ’ 


r-p+ p+ r-p-\- p + 

Rf =-= -jJ^ = a . 


1 - 7i+fi7_i 


the amplitndes of transition and reflection concerning the whole barrier: 


( 12 ) 


At^ E 4 "' = 4 . An= E A^S' = 4-1 ( 13 ) 

n=l n=l 

and the corresponding coefficients of penetrability TmiM and reflection Rmir- 
Tmir = - 7 — Rmir = (14) 

ki 

We check the property 

— \At\^ + \Aii\^ = 1 or Tmir + Rmir = 1, (15) 

ki 

which indicates whether the MIR method gives the proper solntions for the 
wave fnnction. If the energy of the particle is located below the height of step 
with nnmber m, then for a description of the transition of the particle throngh 
the barrier with its tnnneling, the following change mnst be performed: 


km t i ^m) (15) 

where — E). For two rectangnlar barrier steps of arbitrary 

sizes we obtain coincident amplitndes of the wave fnnction calcnlated by MIR 
approach, with corresponding amplitndes obtained in the standard approach 
of qnantum mechanics (np to the hrst 15 digits). Increasing the nnmber 
of steps np into the thousands maintains such an accuracy in coincidence 
between calculations in the MIR approach and the standard approach of 


^We shall analyze the penetrabilities and cross-sections obtained by the multiple inter¬ 
nal reflections method and WKB method. To distinguish the calculated results, we shall 
add abbreviation “MIR” to the “multiple internal reflections method” and WKB to the 
WKB approach. 
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quantum mechanics, and condition fllSp is fulfilled (see Appendix A where 
we briefly present algorithms for the standard approach of quantum mechan¬ 
ics). This is an important test confirming the reliability of the MIR method. 
We have therefore obtained the full coincidence between the solutions for 
all amplitudes obtained by the MIR approach and the standard approach of 
quantum mechanics. A comparative analysis of the accuracy of this method 
and another method in determining the penetrability and reflection coeffi¬ 
cients on the basis of wave functions obtained via the direct integration of the 
radial Schrodinger equation with high accuracy using the Numerov technique 
is given in Appendix E 


2.3. a-nucleus potential and minimization method 

We define interactions between the a-particle and nucleus by the potential 


V{r,l,Q) = vc{r) +UAr(r,Q) +uz(r), 


(17) 


where uc(r), vn{t,Q) and u/(r) are the Coulomb, nuclear and centrifugal 
components. In the spherically symmetric consideration of the a-capture, 
we have 


vc{r) = 


2Ze^ 

r 

Ze^ 


3 - 


at r > r^. 


at r < 


VN{,r,Q) = 


V{A,Z,Q) 


1 -|- exp 


r - m 


and 


= 1.5268+ R. 


(18) 

(19) 


We have previously used parametrization [36|, with a special orientation 
toward the description of a-capture. V (A, Z, Q) is the strength of the nuclear 
component, A and Z are the numbers of nucleons and protons in the nucleus, 
Q is the Q-value of the a-capture, rm is the effective radius of the nuclear 
component, R is the radius of the nucleus, and d is the diffuseness {R, r^, d 
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are in fm). We define the following fnnctions of errors: 


N cr('^"P)(EAr) ’ 

( 20 ) 


where (i?fc) and are theoretical and experimental cross- 

sections of captnre at energy Ek, and the snmmation is performed over ex¬ 
perimental data. We shall search the fusion probabilities po ... Pmax for the 
minimal characteristic fl20|) . We shall call such an approach the minimization 
method. 


£2 = 


£3 = 


a 


1 ^ 

Nti 

1 N |ln(crdh*^°'')(Efc))-ln(a(®^P)(Efc)) 

iV ^ 


(-P)(Efc) 


k=l 


\n{a(^-p){Ek)) 



3. Analysis 

3.1. Dependence of the penetrability on the spatial localization of capture and 
orbital momenta 

For the analysis of the method and general properties of a-capture, we 
shall choose the ^"‘Ca nucleus. Let us consider how the position of rcapture 
influences the barrier penetrability (as we assume capture occurs when the 
maximum of the wave packet describing the a-particle is located within region 
M). The results of such calculations for the energy 2 MeV of the a-particle 
are presented in Fig. [H One can see that small variations in the position of 
^"capture lead to an essential change of the calculated penetrability (for exam¬ 
ple, the influence of the position of rcapture is essentially higher than that of 
the deformation of the nucleus; therefore, there is no sense in studying the 
role of nuclear deformations without a previous consideration of this influ¬ 
ence). The importance of the proper choice of this parameter in calculations 
of the capture cross-sections is now becoming clear. This dependence has an 
oscillating behavior that is characterized by the manifestation of the wave 
nature (origin) of the tunneling processes, and it was previously observed by 
us in the scenarios of proton and a-decays 

As a possible value for rcapture, we choose the coordinate corresponding 
to the minimal amplitude variation of the penetrability oscillations, which 


10 










Figure 1: The penetrability of the barrier depends on the position of the capture 
localization of the a-particle rcapture for the reaction of a Ca at I = 0 and the 
energy of the a-particle = 2 MeV (parameters of calculations: 10000 intervals 
at Tmax = 70 fm, parametrization 


corresponds to the highest stability of the penetrability for possible pertur¬ 
bations (we do not consider similar behavior of the penetrability at r —)■ 0). 
Our performed analysis has shown that such a point is the same for different 
energies of the a-particle, and it corresponds to the minimum of the inter¬ 
nal potential well (we obtain rcapture = 0.44 fm). We shall choose such a 
dehnition of rcapture for further calculations. 

Let us analyze how the penetrability of the barrier depends on the orbital 
momentum. The results of such calculations are shown in Fig. [2] (a). One 
can see that the line of the penetrability is shifted to the right with increasing 
orbital momentum /. This change is explained by increasing the internal well 
of the potential before the barrier, resulting in a larger centrifugal component 
of the potential at larger I (see Fig. [2] (b)). Each subsequent centrifugal 
component of the a-capture cross-section thus gives its own contribution to 
the full spectrum, starting from some larger energy E^. From here, one 
can conclude that the full a-capture cross-section at some sufficiently small 
energy Ea is described by a partial cross-sections at small I only. 

3.2. The fully quantum approach versus the WKB approach in the calculation 
of the cross-sections 

Now, we shall consider how the capture cross-section is determined in the 
fully quantum and WKB approaches. The results of such calculations in the 
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Kinetic energy of a-particle, (MeV) 



Figure 2: The barrier penetrabilities depend on the energy of the a-particle (a) 
and the corresponding potential barriers (b) at different orbital moments I for the 
reaction of a Ca (parameters of calculations: 10000 intervals at rmax = 70 fm, 
parametrization from j^). 



Kinetic energy of a-particle, (MeV) 


Kinetic energy of a-particle, E (MeV) 


Figure 3: The cross-sections of the capture of the a-particle by the ^^Ca nucleus 
calculated in the MIR approach (parameters of calculations: 10000 intervals at 
r ma x = 70 fm, Tcapture Is chosen at the coordinate of the minimum of the potential 
well, which changes depending on the used orbital momentum I, parametrization 
from [i^l. (a) Spectra in the logarithmic scale: the role of contributions at small 
I is more visible, (b) Spectra in the linear scale: the behavior of the summarized 
spectrum after the inclusion of contributions at large I is more visible. One can 
see that each partial cross-section has its maximum. Inclusion of the subsequent 
contributions (at large 1) smooths the full spectrum, slowly transforming it into a 
monotonically increasing curve inside the whole energy range. 
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MIR approach at different values of Rmax are presented in Fig. [3] (here, Fmax 
is the upper limit of the summation in formula ([1]), which we choose instead 
of infinity). The presentation of the cross-sections in the linear scale clearly 
shows their slow decrease after the maxima. Additionally, presenting it on a 
logarithmic scale shows a clear monotonic increase in the cross-sections up to 
the peaks (especially at small 1), which clearly demonstrates the behavior of 
the partial cross-sections at small I on the background of the cross-sections at 
large 1. The WKB approach gives a different result, shown in Fig. 01 We note 



Kinetic energy of a-particie, (MeV) 



Figure 4: The cross-sections of the capture of the a-particle by the “^^Ca nucleus 
calculated in the WKB approach (parameters of calculations: 10000 intervals at 
r ma x = 70 fm, parametrization from [^). In figure (b), one can see the presence 
of sharp peaks in the shape of the spectrum that are hidden on a logarithmic 
scale (see figure (a)). Such a peculiarity of the spectrum is natural for the WKB 
calculations. 


the following two aspects (visible in the linear consideration but hidden in the 
logarithmic): (1) sharp peaks appear after the inclusion of contributions at 
large /; (2) the summarized cross-section is sufficiently higher in comparison 
with the fully quantum calculations presented in Fig. [3l 


3.3. Determination of the probabilities of fusion 

The cross-sections for a-b^'^Ca obtained by the MIR and WKB approaches 
and experimental data 53| are included to Fig. |5l The MIR calculations with 
the inclusion of the fusion probabilities found by the minimization method 
are added to this figure. One can see that the WKB approach gives higher 
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Kinetic energy of a-particle, (MeV) 

Figure 5: The capture cross-sections of the a-particle by the ^'^Ca nucleus obtained 
by the MIR method and WKB approach (parameters of calculations: 10000 inter¬ 
vals at Tmax = 70 fm, parametrization from [s^). Here, the data labeled 1 are the 
experimental data extracted from [s^, dashed blue line 2 is the cross-section at 
^max = 0 short dotted red line 3 is the cross-section at / max = 1, short dash-dotted 
purple line 4 is the cross-section at / max = 5, dash double-dotted orange line 5 
is the cross-section at / max = 10, dashed dark blue line 6 is the cross-section at 
/ max = 12, dash dotted green line 7 is the cross-section at / max = 15, solid brown 
line 8 is the cross-section at /max = 20, dashed dark yellow line 9 is the normalized 
cross-section at /max = 17, solid blue line 10 is the cross-section at / max = 17, and 
dashed dark cyan line 11 is the cross-section at / max = 17 (cross-sections are de¬ 
fined by ([1]) where /max is the upper limit in the summation). The penetrabilities 
are calculated by the MIR method for lines 2-7 and 9-11, and by the WKB ap¬ 
proach for line 8. Lines 10-11 are obtained with the included fusion probabilities, 
and lines 2-9 are obtained without the fusion probabilities. One can see that line 
10, obtained after inclusion of the fusion probabilities, describes the experimental 
data with good accuracy (which is not possible without such coefficients but with 
any arbitrary fitting of the parameter ^|^the a-nucleus potential, see Sect. I3.4jl . 
For line 11, the fusion probabilities are obtained by formulas ()21I) - (I27I1 . 












values for the cross-section in comparison with the calculations from the MIR 
approach at = 10 — 15. We conclude the following: (1) the WKB ap¬ 
proach gives reduced estimations for the fusion probabilities in comparison 
with the MIR approach, and (2) the cross-section obtained in the WKB ap¬ 
proach has discontinuities at higher energies, while the MIR approach gives 
the continuous spectrum shape. 

In the effect of anomalous large-angle scattering (ALAS) was dis¬ 
cussed, which was explained by the sharp angular momentum cut-off ap¬ 
proach at some critical value of (see eqs.(l)-(2) in that paper). This 
case corresponds to calculations by the MIR approach at different values of 
-f^max (where all fusion probabilities are equal to unity at I < Tmax)- Here, 
curve 5 for Lmax = 10 in Fig. [5] (and curve 6 for Lmax = 12) is clearly bet¬ 
ter at describing the experimental data than curve 7 for Lmax = 15 (which 
nearly coincides with subsequent calculations for higher Lmax)- However, all 
these curves are far from the experimental data in comparison with curve 10, 
which includes the fusion probabilities obtained by the minimization method. 
Thus, the sharp angular momentum cutoff approach does not well describe 
the exiting experimental data for a +‘^‘^ Ca, while curve 10 is more success¬ 
ful. This result indicates that the dependence of the fusion probabilities on 
the angular momentum is more complicated and requires more careful study. 
The fusion probabilities used in the calculation of curve 10 are presented in 
Tabl. iBlOl in Appendix B Here, we hnd that po • • • P 4 cire extremely small (in 
contrast to the sharp angular momentum cutoff approach j^). Addition¬ 
ally, pi 5 and pi6 are close to unity (at small pis), which allows one to describe 
small oscillations of the spectrum at high energies (at = 23 — 28 MeV). In 
other words, the presence of the oscillating behavior in the experimental data 
indicates sufficient high fusion probabilities at maximal angular momenta (in 
contrast to the sharp angular momentum cutoff approach [^). 

Eberhard et ah [SSj studied fusion (absorption) of the a-particle by the 
"‘^Ca and '^^Ca nuclei and compared the dependencies of the capture cross- 
sections on the angular momentum 1. These dependencies were obtained 
at one chosen energy (in Fig. 2 in that paper, data are presented for E = 
25 MeV). We shall therefore compare the fusion probabilities for these two 
nuclei as a function of the angular momentum, and such calculations are 
presented in Fig. [6l These data were obtained by analyzing the experimental 
data inside the whole energy region, which indicates a more accurate analysis 
in our approach. Fig. |6] shows that both curves have similar behaviors at low 
I, but they are different at large 1. The fusion probability for '^^Ca is smaller 
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Figure 6: The fusion probabilities for the capture of the a-particle by the nuclei 
^^Ca (solid red line) and ‘^^Ca (dashed blue line) obtained by the MIR approach 
(parameters of calculations: 10000 intervals at r^ax = 70 fm, parametrization 
from 


than that for ^"^Ca at large I - this corresponds to smaller captnre cross- 
sections for ^°Ca than for ^^Ca in the experimental data Decreasing 

the captnre cross-section for ^°Ca at high energies can thns be explained 
by a smaller fusion probability at a large I (these two nuclei have different 
barriers, which give Tmax = 15 for '^^Ca and = 16 for '^^Ca, forming 
an additional contribution to the full spectrum for ^"^Ca). The normalization 
of the calculated spectra on the experimental data (without the inclusion of 
the fusion probabilities) does a worse job of describing the experimental data 
(see curve 9 in Fig. [7j) than curve 10. 

3.4- Optimal parameters of the nuclear component of the potential 

The determination of the potential parameters is a standard task in theory 
of optimization. The analysis shows the presence of only one stable minimum 
in independencies of the function of error on the and d parameters (see 
Fig. [7]), which conhrms the applicability of the optimization methods for 
this task. The determination of these parameters by such a method is given 
Appendix C| in more details. 


m 


In Tabl. [T]we present calculations by the MIR approach with the included 
fusion probabilities and the new parameters found by the optimization tech- 
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Function of error, (mb) 




radius, r [fm] 


Figure 7: Dependencies of the function of errors ei defined by formula (I20p on 
the radius (a) and diffuseness d (b) for the capture a Ca (parameters 
of calculations: 10000 intervals at rmax = 70 fm, applied_MIR method without 
inclusion of the fusion probabilities, parametrization from 


nique for the capture of a Ca and a Ca and for the WKB approach 
without the fusion probabilities for parametrization (such a situation 
did exist before our paper). One can see that with the same parametrization 
used in 3q, the MIR method decreases the error e\ by 1.79 times for a+^°Ca 
and in 2.16 times for a Ca in comparison with the WKB calculations. At 
the new parametrizations found after the inclusion of the fusion probabili¬ 
ties calculations, the MIR method decreases the error e\ by 41.72 times for 
a Ca and by 34.06 times for a Ca in comparison with the WKB cal¬ 


culations at parametrization [3g. An analysis of the Hill-Wheeler approach 
and Wong’s formula shows that they are constructed on an extremely strong 
reduction of the original potential barrier, and even the WKB approach used 
a more correct shape of the barrier in the determination of the penetrability 
(see Appendix D for details). This point is crucial when we extract informa¬ 
tion about the fusion probabilities from the experimental data. The updated 
fusion probabilities for the new parametrization are given in Tabl. [2l These 
results conhrm that the MIR approach (1) is the most accurate at describing 
the tunneling processes through the original barriers and for calculations of 
the fusion probabilities in the a-capture task, and (2) allows one to describe 
the experimental data with the best agreement. 


17 

















a +4^ Ca 


a +4^ Ca 


Parametrization 

own 

Ref. 


own 

Ref. 


Fusion probabilities 

included 

excluded 

included 

excluded 

Penetrability app. 

MIRM 

WKB 

MIRM 

MIRM 

WKB 

MIRM 

1 / 0 , MeV 

28.6444 

28.8374 

36.1 

35.2657 

Tm, fm 

6.5164 

6.6338 

7.3 

6.68670 

d, fm 

0.5107 

0.49290 

0.4349 

0.49290 


0.02141 

0.9233 

0.5154 

0.0246 

0.8401 

0.3881 

£2 

0.02716 

0.5251 

0.3873 

0.0344 

0.4936 

0.3238 

^3 

0.00409 

0.1030 

0.0687 

0.0051 

0.0957 

0.0561 


Table 1: Parameters of the a-nucleus potential and functions of error for the capture 
of a Ca and a Ca. The new parameter determination results for calculations 
of the penetrabilities by the MIR method with the included fusion probabilities given 
in Tabl. m are presented in columns 2 and 5, the MIR calculations of penetrabilities for 
parametrization without the included fusion probabilities are given in columns 3 and 
6, and calculations of the penetrabilities by the WKB approach at parametrization 
without the included fusion probabilities are given in columns 4 and 7. ei, £2 and £3 are 
functions of errors defined by the formula (l20t . 


Nucleus 

Po 

Pi 

P2 

P3 

Pi 

P5 

Pa 

P7 

Ps 

P9 

^^Ca 

0.05 

0.01 

0.01 

0.29 

0.89 

0.94 

0.60 

0.43 

0.83 

0.60 

^^Ca 

0.04 

0.03 

0.01 

0.01 

0.43 

0.31 

0.62 

0.31 

0.71 

1.00 


Nucleus 

Pio 

Pii 

P 12 

Pl 3 

Pii 

Pl5 

P16 

Pn 

P 18 

^'JCa 

1.00 

1.00 

0.01 

0.49 

0.95 

0.01 

0.01 

0.01 

0.01 

Ca 

1.00 

1.00 

1.00 

0.67 

0.09 

0.01 

0.91 

1.00 

0.01 


Table 2: The fusion probabilities for the capture of a Ca and a Ca obtained by 
the MIR approach at parametrization given in Tabl. [1] (see columns 2 and 5 in that table). 
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3.5. The predicted fusion probabilities 

From Fig. E] one can find two aspects that suppress fusion during a- 
capture. 

1) The first aspect is observed at the initial angular momenta. It es¬ 
sentially suppresses the fusion process after the a-particle is transferred the 
barrier and is put into the spatial internal region of the nucleus. This action 
could be nuclear attracting forces between the nucleons forming the main 
contribution of the binding energy of the nucleus and leading to its stability. 
However, this action could also be new forces of dissipative nature with the 
strongest influence on the internuclear motions of nucleons at the smaller 
angular momenta. All such forces should have a unified general form for the 
different nuclei. From Fig. O we find the following dependence describing 
such an aspect: 

^ 1 + eiL-c2)lc3 ’ (21) 

Cl = 1 , C 2 = 4.2, C 3 = 0.5. (22) 

2) The second aspect is shown at the final angular momenta. Comparing 
the curves for the '^°Ca and '^^Ca nuclei in Fig. | 6 l one can find that this 
aspect plays a different role for the two nuclei. Therefore, this aspect should 
be connected with the structure of the nucleus, which can be explained on 
the basis of the closure of nuclear shells. In particular, this aspect can be the 
reason for the difference between the spectra for the '^‘^Ca and ^'^Ca nuclei 
(while the first aspect has the same influence on both nuclei). From Fig. [ 6 l 
we suppose the following form of such an aspect: 

{L-n-Af 

P2{L) = h{L)-Y.e~ C4„ . (23) 

n=l 

Here, the function / 2 (T) suppresses the oscillating dependence with decreas¬ 
ing L, and it can be described as 

/ 2 (L) = (24) 

For the first calculation, we use C 4 „ = 1, C 5 = 0.25 and Cg = 2.5. The param¬ 
eter A should define the influence of the neutrons shells on the calculated 
fusion probabilities, which can be described as 

A = a ■ {N — Amagic) + b. (25) 
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Here, N is number of neutrons of the studied nucleus, and A^magic is fhe 
closest magic neutron number, where A'magic < N (i. e. A^magic = 20 for 
^°Ca and ^"^Ca). We introduce linear interpolation fl25|l on the basis of the 
fusion probabilities given in Fig. |6l Extracting the periods of oscillations 
from Fig. [6l we obtain: 

a = 2.31, b = 4.05. (26) 

The complete fusion probability should be written as 

Pfuii(L) = 1-pi(L)-p2(A). (27) 


However, we recommend such formula for nuclei at close proton shells. 

The fusion probability calculation results for the two studied nuclei “^^Ca 
and ^^Ca by these formulas are presented in Fig. [8l One can see that these 
formulas allow for the sufficient description of the extracted fusion proba¬ 
bilities for both nuclei. The resulting calculation of the cross-section with 
the predicted fusion probabilities for the capture reaction of a Ca is 
added in Fig. EJas the violet dashed line 11. One can see that the calculated 
cross-section with the fusion probabilities, when predicted in such a manner, 
describes the experimental data better than any other calculation without 
consideration of fusion during the a-capture. 

The ^°Ca nucleus is double magic, while the "^^Ca nucleus only has a 
magic number for protons Z = 20. Therefore, the binding energy for the 
^°Ca nucleus is higher than that for “^^Ca. Therefore, the ^°Ca nucleus is 
more stable and it is more difficult to synthesize a new one by means of 
the scattering of the a-particle off this nucleus. In particular, the fusion 
between the a-particle and this nucleus during such a scattering reaction 
should be maximally weakened and, in some cases (for some orbital momenta, 
etc.), even vanishing. This situation is clearly demonstrated by our results 
presented in Fig. [6l where one can see smaller fusion probabilities for ^°Ca 
than for ^^Ca at higher angular momenta. This situation is also conhrmed 
by the experimental cross-sections for a Ca and a Ca. 


In 1975, Thibault et al [73 


observed the simultaneous coexistence of two 
types of shapes in the magic ^^Na nucleus - spherical and deformed in the 
ground state. The presence of the deformed shape in the ground state was 
explained by the strong correlation between 2p-2n excited energies for the sd 
and p/-shells, which results in an increase of the binding energy and stability 
reinforcement for the nuclei near the shell at A^ = 20. In subsequent studies. 
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Figure 8: The predicted fusion probabilities for the capture reactions of a Ca 
(a) and a+'^^Ca (b). Here, solid blue lines are calculations for the predicted fusion 
probabilities given by formulas (I2ni - (l271) . and dashed red lines are the extracted 
data of the fusion probabilities given in Fig. [H 

it was found that such a coexistence of two shapes is observed for majority 
of the nuclei near the shells at = 20 and N = 28 j^. 

In this connection, we suppose that each shape (in the ground state) 
determines the most stable own state of the nucleus, and it should therefore 
have a definite relation with the suppression of the fusion processes if one 
considers the scattering of the a-particles on such nuclei. In particular, both 
shapes should appear in the nuclei with a close neutron shell at iV = 20, 
which one can look for in the experimental data of the a-capture reactions. 
^°Ca is such a nucleus, for which we see the more essential occurrence of 
the second suppression aspect of fusion at higher energies of the incident 
a-particles (see Fig. [6]). This process leads to a more essential weakening of 
the a-capture cross-section for such a nucleus at higher energies. 

Each shape should be related to its own type of fusion suppression in 
the a-capture reactions for the ^°Ca and '^^Ca nuclei. In particular, the 
spherical shape of the nucleus does not involve higher angular momenta, 
so it should be related to the hrst aspect (shown at the smallest angular 
momenta). The deformed shape can be related to higher angular momenta 
(that corresponds to higher energies in cross-sections), which is connected 
with the second aspect of reducing the fusion. We observe the coexistence 
of two such peculiarities in our results in Fig. |6l which are correlated with 
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Fusion probability, 


the experimental data for both nuclei. By such logic, a deeper study of the 
coexistence of the spherical and deformed shapes of the nuclei in the ground 
state can be recommended via reactions of the a-capture by the interesting 
nuclei (both theoretically and experimentally). 

Conhrmation of the coexistence of the two shapes for the nuclei near 
the shells at iV = 20 and N = 28 has lead to the reconsideration of the 
main positions of the nuclear shell model and to the discovery of new magic 
numbers A review of this topic ^ gives interesting indications about 


the new neutron magic numbers at = 16 and N = 26 and the properties 
of such nuclei (while the standard theory gives us only seven experimentally 
known neutron numbers at 2, 8, 20, 28, 50, 82, 126). In this regard, it could 
be interesting for experimentalists to propose to investigate the fusion process 
at the capture of the a-particle by the ^®Ca nucleus. Such information could 
provide new insight into our understanding of physics of nuclei with such a 
neutron magic shell. We support such research by our predictions for this 
nucleus given in Fig. |9l 



4 6 8 10 12 14 16 18 20 22 24 26 28 30 



Kinetic energy of a-particle, E (MeV) 


Figure 9: The fusion probabilities calculated by formulas (l2Tl) - ([?7l) (a) and cross- 
sections (b) for capture of the a-particle by the ^°Ca, ^‘^Ca and “^^Ca nuclei obtained 
by the MIR method (parameters of calculations; 10000 intervals at r^ax = 70 fm). 
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4. Conclusions 


In this paper, the multiple internal reflections method is generalized in 
the description of the capture of a-particles by nuclei. On this basis, the 
a-captures by the ^°Ca and ^"‘Ca nuclei are analyzed. Note the following: 

1. Comparing the capture cross-sections calculated by the MIR method 

we And (1) new fusion probabilities 


with the experimental data 


of the a-particles with the “^^Ca and '^^Ca nuclei (see Tabl. |2]), and (2) 
new parameters of the a-nucleus potential (see Tabl. [T]). Considering 
the fusion processes and including them in the calculations essentially 
allows us to increase agreement with experimental data (see Tabl. [T]). 
In particular, the MIR method found new parametrizations with the 
included fusion probabilities and decreased the error Si by 41.72 times 



nations). A more realistic picture is given in Table lBl3] (corresponding 
to curve 10 in Fig. [7]). In particular, the fusion probabilities at the first 
angular momenta are close to zero, contradicting the idea of a sharp 
angular momentum cutoff (see Tabl. [2]). 

3. It is shown that (1) the WKB approximation gives reduced probabili¬ 
ties of fusion in comparison with the MIR approach, and (2) the cap¬ 
ture cross-section obtained in the WKB approach has discontinuities 
at higher energies (the MIR approach provides a continuous spectrum 
shape). 

4. It is shown that Wong’s formula and the Hill-Wheeler approach for 
barrier penetrability determination (1) use essentially a reduced inverse 
parabolic approximated shape of the original barrier, (2) do not use 
the shape of the realistic barrier (with the exception of the localized 
region near the barrier maximum point), and (3) do not correspond to 
standard penetrability in quantum mechanics. 

5. A new formula of the fusion probabilities is proposed. Explanations of 
the differences between the cross-sections for a Ca and a 4-^'^ Ca 
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and the previously discovered coexistence of the spherical and deformed 
shapes in the ground state for nuclei near the neutron magic shell N = 
20 are given by our formula for fusion probabilities. To obtain deeper 
insight into the physics of nuclei with the new magic number N = 
26, a cross-section of a Ca is predicted for further experimental 
conhrmation. 

6 . We compared the effectiveness of the MIR approach in determining 
of the penetrability and reflection coefficients with an approach based 
on the determination of the wave functions obtained via the direct 
integration of the radial Schrodinger equation with high accuracy using 
the Numerov technique. For analysis, we chose the reaction a Ca 
for case / = 0. We conclude the following (see Appendix E): 

(a) To estimate the accuracy of the determination of the penetrabil¬ 
ity and reflection coefficients, we use the standard test of quan¬ 
tum mechanics, calculating parameter £i = | |T -|- i?| — 1| by each 
method (here, T and R are the penetrability and reflection coeffi¬ 
cients). The Numerov approach determines these coefficients with 
accuracy Si = 10“^, while the accuracy of the MIR approach is 
restricted by computer capability (see Tables IE. 71 and IE.81 where 
we provide the MIR calculations at ei < 10“^® for any chosen r^ax 
and N). 

(b) Analyzing the stability of determination of the penetrability and 
reflection coefficients after varying the parameters r^ax and N, we 
conclude that the MIR approach is essentially more accurate than 
the Numerov approach. The accuracy of determination of these 
coefficients by the Numerov approach is restricted by limit 10“^: 
we have only 3 stable digits for these coefficients at variation of 
Tmax and number of intervals N. Simultaneously, the accuracy of 
the determination of these coefficients by the MIR approach is 
inhnitely increased by increasing the number of intervals N: we 
reach the hrst 8 stable digits of these coefficients when N —)■ 50000 
(see Tables lETl and [ETSjl . 

In summary, we conclude that the MIR approach is more accurate and 
stable. Simultaneously, it requires a similar time for calculations as the 
Numerov approach. 

This result can be explained by the following. The MIR approach gives 
exact analytical solutions for the wave function (for the chosen bound¬ 
ary conditions) for the potential with any hnite number of rectangular 


Appendix E 


24 











steps. The accuracy of the description of the original potential by the 
rectangular steps potential is determined only by the number of steps, 
where the limit of increasing the number of such steps is restricted only 
by computer capability (i.e., in principle, the methodological accuracy 
of such an approach can be increased up to inhnity by the presented 
formalism in this paper). 

At the same time, the accuracy of the Numerov approach has some 
limits that cannot be overcame. The accuracy of such an approach is 
restricted by the following: (1) The Numerov technique gives values 
for the wave function in the next radial coordinate on the basis of the 
values of the wave function in the previous two radial coordinates with 
some error, i.e., it is approximated. Such accuracy is only valid inside 
a sufficiently small radial region. (2) Error in the determination of the 
wave function inside the whole radial region of action of nuclear forces 
on the basis of the Numerov technique is increased (with increasing 
number of intervals N). (3) The amplitudes of the Coulomb functions 
at the near asymptotic region are calculated on the basis of asymptotic 
series that are divergent (usually, such functions are used at the best 
convergent limit). 

7. We analyzed the role of the nuclear deformations in the determinations 
of the fusion probabilities in the a-capture task for a Ca. We hnd 
the following (see Appendix F): 

(a) Analysis of changes in the penetrability of the potential barriers 
after the inclusion of the deformation of the nucleus for case I = 0 
provides the clearest understanding of such a question. The largest 
influence of the nuclear deformation exists for angles of deforma¬ 
tion 6 of 0° and 90° and for the lowest energy of the a-particle, 
while at 6 = 55°, such an influence is missing (see Fig. IF.12j) . 

(b) The MIR approach allows one to obtain the fusion probabilities 
for the deformed nuclei and to describe the experimental data of 
the a-capture cross-sections with high accuracy similar to that for 
the spherical nuclei (see Tabl. IF.91 and lines 2 and 3 in Fig. IF.13p . 
We observe the presence of small oscillations in the calculated 
spectrum after the inclusion of the nuclear deformation. 

Another important aspect is the determination of the boundary where 
the transmitted flux is defined. It turns out that small displacements of 
this boundary essentially change the resulting penetrability and reflec¬ 
tion coefficients. For example, for a Ca at / = 0 and = 2 MeV 
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the role of such a displacement can be estimated by the ratio fi = 
^ (^capture 7.42 fm)/T(reapture = 4.5 fm) = 4.09869 ■ 10-V2.07705 ■ 
10“^ = 1.9733 (see Fig. 1). Simultaneously, the role of nuclear de¬ 
formation in this task can be estimated by the ratio of the pene¬ 
trabilities for spherical and deformed nuclei for this reaction. We 
have (Tdef — Fspher)/Fspher = 0.58015 for Ea = 2.1 MeV and I = 0 
(see Fig. F.12 (b) at close energies). From here, we obtain the ratio 
/2 = Fdef/Tspher = (Fdef — Fgpher)/Fgpher + 1 = 0.58015 -|- 1 = 1.58015. 
One can see that this ratio /2 is less than the ratio /i. From here, one 
can conclude that in fully quantum calculations, the role of the position 
of the boundary Tcapture in determining the resulting penetrability (and 
cross-sections) is more important than the nuclear deformations. To 
study such an aspect correctly in the a-capture task (and to work with 
it), some formalism should be constructed. We use the MIR approach 
for this aim. This is one of reasons to use the MIR technique. 

We summarize that the fusion processes are an essential ingredient in nuclear 
capture, and we propose the use of the MIR method as a fully quantum highly 
precision tool for the study of such processes in nuclear tasks. 
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Appendix A. Direct method 

We shall soon add a solution for the amplitudes of the wave function 
obtained by the standard technique of quantum mechanics that could be 
obtained if used only in the condition of continuity of the wave function 
and its derivative at each boundary, but on the whole region of the studied 
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potential (key relations were given in such an approach was further 
applied for the study of a-decay, proton decay [^, and tunneling in quantum 
cosmology tasks [70|, ITIJ). First, we hnd the functions fN-i and gN-i (from 
the last boundary at tn-i)- 


_ kN-l + kN 

jN-i — -j e , Qn-i — -j , 

kN-i — k^ k^ — k]si-i 


N gi(A:iv+fciv-i)»'iv-i 


Then, using the following recurrent relations: 


fj-i = 


{kj_i - kj) + fj {kj_i + kj) 


2ikj—irj — i 


(A.l) 


(A. 2 ) 


we calculate the subsequent functions /v-2, /at-s, /v-4 ■ ■ ■ fi, and by such 
a formula: 


9j-i — 9j 


2kj QKkj-l+kj)rj-i 


(A.3) 


(fcj_i + kj) + fj {kj_i - kj) 

the functions gN- 2 , gN-s, fi'Af -4 ■ ■ ■ gi- From fi and gi, we hnd amplitudes 
Q?! and /Si and the amplitude of transmission At'- 


Q?! — 0 , At — /Si — 


9i 

'fi 


(A.4) 


Now, using the recurrent relations: 


f3j+i 


(3j + aj 


9j+i ^ 




(A.5) 


and such a formula: 

Oij = fj+9j, (A. 6 ) 

we consistently calculate the amplitudes 0:2, /S 2 , 0:3, /S 3 ...aN-i, /Sat-i- Fi¬ 
nally, we hnd the amplitude of rehection Aji: 


Ar = /Sjv-i + aN-i _ fnkNTN-i 


(A.7) 


As a test, we use condition: 

7 — |AtP + \^r\^ = 1- (A. 8 ) 

kn 

Studying the problem of proton decay, we used a technique to check the am¬ 
plitudes obtained previously by the MIR approach and obtained coincidence 
up to the hrst 15 digits for all considered amplitudes. The results for the de¬ 
pendence of the penetrability of the position of point r^apture and the obtained 
cross-sections are thus independent of the fully quantum method used. 
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Appendix B. Calculations of the fusion probabilities 


In Tabl. IB.31 we present the resnlts of the fusion probabilities calculations 
for the capture reaction a+"^^Ca, with in dependence on the different maximal 
values of the orbital momentum Lmax (and the corresponding functions of 
errors). Here, the penetrabilities of the barrier are calculated by the MIR 


method, and the parametrization of the a-nucleus potential is used from |36 
From such results, one can see that at a smaller maximal momentum Lmax all 
probabilities are equal to unity. However, at increasing Lmax starting from 
-hmax = 9 5 the fusion probabilities at the hrst angular momentum values 
are close to zero. This tendency of the vicinity to zero of the hrst fusion 
probabilities is present at increasing Lmax up to the maximum of its possible 
value. Such a tendency is also present for another reaction a Ca, after 
hnding its own parametrization of the a-nucleus potential (see Tabls. [1] and 
12]). Such a picture clearly demonstrates the failure of the idea for a sharp 


angular momentum cutoff approach introduced by Glas and Mosel [63|, l64 
A more accurate description of the tendencies of the fusion probabilities 
at smaller and larger angular momenta allows us to describe the difference 
between the experimental data 53| for reactions a Ca and a Ca. 


This indicates the importance of cancelation of the sharp angular momentum 
cutoff approach (which has produced countless papers) in the capture and 
synthesis tasks if we want to study the fusion processes. 


Appendix C. Calculations of the potential parameters by the min¬ 
imization method 

To estimate errors given by the different methods for the capture of 
a Ca and a Ca, such calculations should be compared for the same 


parametrization. We choose the parametrization from (3^, as this paper 


was directly oriented to the study of o-capture processes, and the authors 
provide the parameters for such reactions. We include into the analysis our 
MIR approach, the WKB approach and the approaches developed by Hill 
and Wheeler in and Wong in to calculate the penetrability for the 
inverse parabolic barrier, which approximates the original potential barrier 
(we shall call them Hill-Wheeler and Wong’s approaches). From the results 
presented in Tabl. 1C.41 one can see that the MIR method has half as many 
errors as the WKB calculations, while the Hill-Wheeler and Wong’s formulas 
give estimations of errors that are similar to those of the MIR approach. 
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Po 
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0.01 

0.005 

0.156 

Pi 

1.0 

1.0 
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1.0 

0.16 

0 

0 

0 

0.01 

0 

0.001 

0.090 

P2 

1.0 

1.0 

1.0 

1.0 

1.0 

0 

0 

0 

0.02 

0 

0.001 
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P3 
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1.0 

1.0 

0 

0 

0 

0.02 

0.02 
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0.204 

P4 

1.0 

1.0 

1.0 

1.0 

1.0 

0.58 

0 

0 

0.26 

0.19 
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P5 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

0.73 

0.45 

0.86 

0.81 
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P6 

- 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

0.22 

0.36 
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0.455 

Pr 

- 

- 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 

1.0 
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- 
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- 
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1.0 

Pn 

- 
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0.411 

0.347 

0.283 

0.209 

0.129 

0.0569 

0.0364 

0.0314 

0.0301 

0.0293 

£2 

1.962 

1.254 

0.825 

0.545 

0.363 

0.298 

0.172 

0.0735 

0.0472 

0.0420 

0.0405 

0.0397 

£3 

0.171 

0.125 

0.095 

0.071 

0.054 

0.041 

0.024 

0.0107 

0.0070 

0.0062 

0.0060 

0.0059 

■^min 

-16.6 

-14.1 

-11.3 

-8.26 

-4.90 

-1.29 

2.56 

6.50 

10.9 

15.5 

20.24 

25.15 

-^max 

8.45 

9.41 

10.5 

11.8 

13.4 

15.5 

17.1 

19.3 

21.7 

24.3 

27.32 

30.52 


Table B.3: Fusion probabilities calculation results for reaction a Ca for parametriza- 
tion by the minimization method, where the MIR approach is used in calculations of 
the penetrabilities, po ■ • - Pig are the fusion probabilities, ei, £2 and £3 are functions of 
errors defined via the formula ()20p . Flmin is the minimum of the potential well before the 
barrier (MeV), F'max is the height of the potential barrier (MeV), and Lmax is the maximal 
value of the orbital momentum, used in the summation in eq. ©■ 
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Method 

a +4^^ Ca 

a +44 Ca 


£2 

£3 


^2 

£3 

MIR method 

0.5154 

0.3873 

0.0687 

0.3881 

0.3238 

0.0561 

WKB approach 

0.9233 

0.5251 

0.1030 

0.8401 

0.4936 

0.0957 

Hill-Wheeler app. 

1.3923 

0.6236 

0.1306 

0.3927 

0.3265 

0.0569 

Wong formnla 

0.6834 

0.4527 

0.0829 

0.5507 

0.4007 

0.0713 


Table C.4: Errors obtained by different approaches for cross-sections of the capture reac¬ 
tions a Ca and a Ca with the potential parametrization given in [^ . where we 
do not include calculations of the fusion probabilities. Si, £2 and £3 are functions of the 
errors defined by the formula (l20)) . From the presented results, one can see that for the 
same parametrization [s^ , the MIR method decreases errors twice in comparison with the 
WKB calculations, while the Hill-Wheeler and Wong’s formulas give estimations of errors 
that are similar to those of the MIR approach. 


From the first consideration, it strangely was not clear why the difference 
between the resnlts obtained by the MIR approach and Hill-Wheeler and 
Wong’s formnlas was so small. The MIR approach essentially inclndes more 
formalism (and compnter resonrces) for the approximated potential barrier 
and calcnlations of its penetrability than the Hill-Wheeler and Wong’s formn¬ 
las. However, from snch a point, it has become clear that the errors are only 
the compared difference between the calculations and the experimental data 
of the cross-sections, and they absolutely ignore the level of accuracy of the 
approximation of the original barrier. According to snch logic, we can take 
another potential, calcnlate its penetrability by any approximated formnla, 
and, if this new valne is closer to the experimental data, conclnde the snccess 
of the nsed approach. However, this is absolntely the wrong conclnsion. The 
error shonld therefore also take into acconnt the level of the approximation 
of the original barrier by the method nsed. As the WKB approach has worse 
accnracy, we do not inclnde it in any fnrther analysis. 

Let ns see which parameters of the a-nnclens potential the stndied meth¬ 
ods give. The resnlts obtained by the optimization techniqne (we choose the 
Gradient descent method) for a + ^Ala are presented in Tabl. 1C.51 where 
we do not inclnde calcnlations of the fnsion probabilities. One can see that 
the estimations for the radins, rm, are very close for any nsed approach with 
parametrization 3^. The strength of the potential, Vq, obtained by the MIR 
approach and Wong’s formnla are close to parametrization j^, in contrast to 
the Hill-Wheeler approach. All three approaches essentially give different es- 


30 

















Parameters of potential 

errors 


Vo, Mev 

Tm, fm 

d, fm 

-^min 

-^max 

£1 

£2 

£3 

MIRM 

34.091 

6.6855 

0.1789 

-24.04 

7.48 

0.0738 

0.0855 

0.0129 

HWA 

16.463 

6.6866 

0.1568 

-6.67 

7.71 

0.0761 

0.0987 

0.0151 

WoF 

33.155 

6.6864 

0.1455 

-23.30 

7.66 

0.0730 

0.0879 

0.0129 


Table C.5: Parameters of the potential obtained via minimization of the function of the 
error Si for the capture reaction a Ca, where we do not include calculations of the 
fusion probabilities, ei, 62 and £3 are error functions defined by formulas (| 20 )) . illmin is 
the minimum of the potential well before the barrier (MeV), and Emax is the height of 
the potential barrier (MeV). MIRM, HWA, and WoF are the MIR method, Hill-Wheeler 
approach, and Wong’s formula, respectively. One can see that after the minimization pro¬ 
cedure, the errors are less than a quarter of those in the previous calculations in Tabl. 1C.41 
However, at a similar accuracy of experimental data description, different methods give 
essentially differences between parameters Vq and d. From here, it follows that the level 
of agreement with the experimental data on the basis of the given errors is not sufficient 
for extracting reliable information about the parameters of the interacting potential. 


timations for diffuseness, d, in comparison with the results from j3^ . It now 


becomes clear that the included formalism with different approximations of 
the original barrier gave different estimations of the parameters at similarly 
obtained errors. This result underlines the importance of including the level 
of the approximation of the initial barrier when determining the error of the 
method. 


Appendix D. Comparison of the approach of Hill and Wheeler 
and Wong’s formula vs. the MIR and WKB ap¬ 
proaches 

The direct comparison of the calculations from the MIR approach and the 
Hill-Wheeler formula shows that the penetrability obtained by the Hill and 
Wheeler approach is essentially higher at smaller energies (see Fig. ID.lOp . 
but this overestimated result has no visible influence on the cross-section. 
However, this result is not correct in a fully quantum consideration, where 
we essentially have an influence of the shape of the barrier on the estimation 
of the penetrability. To clarify the role of the shape of the inverse parabolic 
barrier (used in the Hill-Wheeler approach), we have calculated the penetra¬ 
bility of such a barrier by the MIR method. The result of such calculations 
is presented in Fig. ID.Ill (see dash-dotted green line in this hgure), which is 


31 





















Penetrability of barrier 



Kinetic energy of a-particle, E (MeV) 



Figure D.IO: The penetrabilities of the barrier for the reaction a Ca at / = 0 
calculated by the MIR approach (solid blue line) and by the Hill-Wheeler formula 
(dashed red line) (parameters of calculations: 1000 intervals at r^ax = 70 fm). On 
the linear scale, we obtain a sufficiently close agreement between the penetrabilities 
calculated by the two approaches (see figure (a)). However, a more accurate study 
of the logarithmic scale shows a higher estimation for the penetrability obtained 
by the Hill-Wheeler formula (see figure (b)). Now, it becomes clear that the 
application of such a formula for the determination of penetrabilities in the nuclear 
decay tasks can lead to essential errors in estimations of half-lives. In fusion tasks, 
the Hill-Wheeler formula is saved by the fact that for the determination of the 
capture cross-section, we need to summarize up to 18 curves of the penetrabilities 
with shifted maxima (only where this formula is working). 
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strongly different from the estimation of the Hill-Wheeler formula (see dashed 
red line in this hgure). Only the penetrability at the point of intersection 



Figure D.ll: The penetrability of the inverse parabolic barrier used in the approach 
of Hill and Wheeler in at Z = 0 for reaction a Ca as calculated by the 
MIR method (parameters of calculations: 1000 intervals at rmax = 70 fm). Such 
calculations by the MIR approach are presented by the dash-dotted green line. 
For comparative analysis, we add calculations of the penetrability of the original 
barrier by the MIR approach (see solid blue line) and by the Hill-Wheeler formula 
(see dashed red line), shown in the previous figure. This result clearly shows that 
the MIR method essentially gives a different penetrability than the Hill-Wheeler 
formula. The reason for such a difference between the calculations obtained by 
the two approaches will become clearer if we compare the shapes of the inverse 
parabolic barrier with that of the original barrier. Such a comparison shows that 
there are no similarities, with the exception of a small neighborhood near the 
barrier maximum (which gives close estimations of the cross-sections). 

between all the curves in such a hgure plays a role in the calculations of the 
cross-section, which explains the not bad agreement of the calculations by 
the approach of Hill and Wheeler with the experimental data. However, from 
the analysis and calculations above, it is clear that good agreement of the 
calculations by the approach of Hill and Wheeler with the experimental data 
does not conhrm the validity of such a formula. In Tabl. ID.6( we demonstrate 
the convergence of calculations of the penetrability by the MIR approach for 
such an inverse parabolic barrier. Moreover, the WKB approach includes the 
shape of the original barrier in the calculations. It is therefore more correct 
and accurate in the estimations of the penetrabilities of the barriers in fusion 
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Number of 

r 

max = 50 fm 

?^max = 200 fm 

intervals 

1 MeV 

15 MeV 

30 MeV 

1 MeV 

15 MeV 

30 MeV 

10000 

0.48634 

0.53097 

0.57703 

0.49473 

0.51760 

0.54100 

20000 

0.48086 

0.52368 

0.56866 

0.44806 

0.49275 

0.54366 

30000 

0.48226 

0.52432 

0.56828 

0.47961 

0.51829 

0.55901 

40000 

0.48018 

0.52204 

0.56617 

0.48752 

0.53010 

0.57372 

50000 

0.48019 

0.52197 

0.56601 

0.47070 

0.51346 

0.55993 

(^max) 

-3747000 

-95000000 


Table D.6: The penetrability of the inverse parabolic barrier for a+'^'^Ca at 1 = 0 as calcu¬ 
lated by the MIR method for different values of the following calculated parameters: The 
number of intervals and external boundary r^ax- The parameters of the inverse parabolic 
barrier are determined according to the Hill-Wheeler approach (6(i |. which approximates 
the original potential barrier with parametrization [s^. In such calculations we do not 
include the fusion probabilities. H(rinax) is the inverse parabolic potential at point rmax 
(MeV). 


tasks (in comparison with Hill-Wheeler approach and Wong’s formula). 


Appendix E. Comparison of the MIR approach vs. the Numerov 
approach 

The penetrability and reflection coefficients can be determined on the 
basis of wave functions obtained via the direct integration of the radial 
Schrodinger equation with high accuracy using the Numerov technique. In 
this section, we shall compare the effectiveness of such a method in deter¬ 
mining the studied coefficients in comparison with the MIR approach. For 
simplicity of analysis, we shall study the reaction a -1-"^^ Ca for case / = 0. 

We use the same separation of the radial region on A regions as dehned 
in (E]), and we assume that the capture of the a-particle with the nucleus 
occurs in the hrst region after its tunneling through the barrier. We dehne 
the normalized radial wave function as 


X{r) 


exp{—ikir), at < r < ri (region 1), 

A 2 C 2 (r)-b i ?2 d 2 (r), at ri < r < 

^ max: 


(E.l) 


where A 2 and B 2 are unknown complex amplitudes, and C 2 (r) and d 2 {r) 
are linear independent partial solutions of the radial Schrodinger equation 
with potential dehned by ffT7)l -ffT ^ . These amplitudes should be found from 
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continuity conditions of the full wave function and its derivative at point ri. 
As a result, we obtain: 


3 f?(A 2 ) = Nq [cos(/ciri) ^2(^1) + kism{kiri) d2{ri)], 
A(A2) = No [ki cos(/ciri) d2{ri) - sin(/ciri) 

3 f?(fi 2 ) = - No [cos(/ciri) ^(ri) + /ci sin(/ciri) C2(ri)], 
^>(^2) = - No [ki cos(/ciri) C2(ri) - sin(A;iri) ^(ri)], 


where 


No = 


1 

C 2 (ri) 4(^^i) - C2(n) d 2 (ri) ‘ 


(E.3) 


We dehne the partial wave functions C 2 (r) and d 2 {r) so that in the asymp¬ 
totic region, starting from some selected boundary they corresponds to 
Coulomb functions as C 2 (r) = Gi=o(r), d 2 (r) = Fi=o(r). To determine the un¬ 
known wave functions and their derivatives in region ri < r < rag (where the 
nuclear component of the potential is important), we calculate them numeri¬ 
cally using the Numerov approach. One can rewrite the full wave function in 
the asymptotic region via a combination of the incident and reflected waves 
relative to the barrier as 


X(r) = Ainc [Go(r) - i Fo(r)] -f- [Go(r) + i Fo(r)], (E.4) 


where we have introduced the amplitudes Ai^c and Aref for the incident wave 
and reflected wave. Such amplitudes can be rewritten as 

^ _ A2 + i B2 ^ _ A2 — i B2 

^inc 2 ’ ^ref ^ 

and we have the amplitude of transition equaling unity. At,. = 1. Taking into 
account the property Go Fq — GqFo = 1 for the Coulomb functions, we obtain 
the following solutions for the penetrability and reflection coefficients: 



4fci 1 D ^ /i - h 
k /1 + /2’ /1 + /2’ 


(E.6) 


where 

ft = [5J(A2)]^ + [A(A2)]2 + mB2)? + [A(i?2)]^ 

/2 = 2 [A(A2) 3fJ(52) - 3fJ(A2) %B2)]. ^ ^ 

To compare the effectiveness of the Numerov and MIR methods in deter¬ 
mining the penetrability and reflection coefficients, we calculate these coeffi¬ 
cients in dependence on the external boundary rmax and number of intervals 
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’^max, fni 

Numerov approach 

MIR a] 

pproach 


T 

R 


T 

R 

50 

0.890661 

0.109292 

4.5 X 10“^ 

0.99997095707 

2.90429 X 10-5 

75 

0.890581 

0.109371 

4.6 X 10-® 

0.99997116245 

2.88375 X 10-5 

100 

0.890661 

0.109292 

4.5 X 10-5 

0.99997139532 

2.86046 X 10-5 

150 

0.890822 

0.109133 

4.4 X 10-5 

0.99997145272 

2.85472 X 10-5 

200 

0.891300 

0.108659 

3.9 X 10-5 

0.99997147408 

2.85259 X 10-5 

250 

0.890982 

0.108974 

4.2 X 10-5 

0.99997142915 

2.85708 X 10-5 


Table E.7: The penetrability T and reflection R of the barrier with parametrization 
for a Ca at / = 0 and energy E = 25 MeV calculated by the Numerov approach and 
MIR method, for different values of the external boundary r^ax- Here, ei = ||r + i?| — 1| 
is the error calculated to test the obtained coefficients; for all MIR calculations, we have 
Si < 10“^®, we do not include the fusion probabilities, and the number of intervals is the 
same and equals 10000. 


Number of 
intervals 

Numerov approach 

MIR a] 

oproach 

T 

R 

Si 

T 

R 

100 

0.934547 

0.065452 

8.9 X 10-^ 

0.95572409894 

0.044275901 

1000 

0.894080 

0.105897 

2.1 X 10-5 

0.99997192337 

2.80766 X 10-5 

10000 

0.890822 

0.109133 

4.4 X 10-5 

0.99997145272 

2.85472 X 10-5 

25000 

0.890726 

0.109228 

4.5 X 10-5 

0.99997144515 

2.85548 X 10-5 

50000 

0.890533 

0.109419 

4.7 X 10-5 

0.99997144338 

2.85566 X 10-5 


Table E.8: The penetrability T and reflection R of the barrier with parametrization 
for a +‘^'^ Ca at Z = 0 and energy E = 25 MeV calculated by the Numerov approach and 
MIR method for different values of the number of intervals. Here, ei = ||T + i?| — I| is 
the error used to test the calculated coefficients, and for all MIR calculations, we obtain 
El < 10“^®. The external boundary is Cmax = 150 MeV, and in such calculations we do 
not include the fusion probabilities. 
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N by both approaches. The results of these calculations are presented in 
Tables IE.71 and IE.81 According to such results, we conclude the following: 

1. To estimate the accuracy of the determination of the penetrability and 
reflection coefficients, we use the standard test of quantum mechanics, 
calculating the parameter Si = ||T + i?| — 1| by each method. The 
Numerov approach determines these coefficients with accuracy ei = 
10“®, while the accuracy of the MIR approach is restricted by computer 
capability (in the Tables above, we provide the MIR calculations at 
El < 10“^® for any chosen rmax and N). 

2 . Analyzing the stability in determining the penetrability and reflection 
coefficients after varying the parameters r^ax and N, we conclude that 
the MIR approach is essentially more accurate than the Numerov ap¬ 
proach. The accuracy of the determination of these coefficients by the 
Numerov approach is restricted by the limit 10“^: we have only 3 stable 
digits for these coefficients at the variation of r^ax and number of inter¬ 
vals N. At the same time, the accuracy of the determination of these 
coefficients by the MIR approach is inhnitely increased by increasing 
the number of intervals N: we reach the hrst 8 stable digits of these 
coefficients when N —>■ 50000 (see Tables lETl and lEl^ . 


Appendix F. Role of the nuclear deformations in the a-capture 
task 


In this section, we shall analyze how the obtained results are changed 
if we take into account the nuclear deformation of the target-nucleus. In 
this paper, we shall consider the case of axial deformation of the nuclei. For 
the ^nucleus interactions, we use the potential with parametrization dehned 
which was found for the o-capture task with deformed nuclei. The 
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cross-section for o-capture by nuclei with axial deformation can be dehned 
as 


cr, 


capture 


{E) = 


71 h 


2 '"i‘^+00 


2m E 


^(2/ + l)Ti(0)Pi sin Ode, 


(F.l) 


1=0 


where 9 is the angle between the axial symmetry axis of the deformed nucleus 
and the line connected centers-of-masses of the a-particle and the nucleus, 
Ti{6) is the penetrability of the barrier, which depending on the angle 9, and 
Pi is the probability for fusion of the a-particle and nucleus. 
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For simplicity of the analysis, we shall study the capture reaction of 
a Ca. The experimental data on the static quadrupole deformation 
parameter /? is taken from the RIPL-2 database [7^ (for ^'^Ca we have 
/3 = 0.175423). First, we shall analyze how the penetrability of the bar¬ 
rier is changed after the inclusion of the deformation of the nucleus. In 
Fig. IF.121 (a), we present our calculations for the penetrability, as a function 
of the angle of deformation 9 for case I = 0. We see that the largest influence 




Angle of deformation of nueleus, $ (degree) Kinetie energy of a-particle, (MeV) 

Figure F.12: The ratio of the penetrabilities of the barriers for the deformed nucleus 
vs. the spherical one, as a function of the angle of deformation 6 at the incident 
energy of the a-particle of 25 MeV (a) and the energy of deformation 6 at the angle 
of deformation of 9 = 90° and (b) for reaction a Ca at / = 0 calculated by 
the MIR method (parameters of calculations: 10000 intervals at rmax = 150 fm). 
Here, Tjef is the penetrability of the barrier for the deformed nucleus, Tspher is the 
penetrability of the barrier for the spherical nucleus. 

of the nuclear deformation exists for angles of 0° and 90°, while at 6* = 55°, 
such an influence is missing. Choosing the angle 6 = 90°, we can then ana¬ 
lyze the role of nuclear deformation in the calculations of the penetrability, 
as a function of the energy of the incident a-particle. The results of such 
calculations are presented in the subsequent Fig. IF.121 (b). From here, one 
can see that high importance of taking into account the nuclear deformations 
for the lowest energy of the a-particle (for = 5.1 MeV, where we obtain 
a ratio of (Tdef - Tspher)/Tspher = 0.27739). 

In Tabl. IF.91 we present the results of the determination of the fusion 
probabilities by the MIR approach for a Ca after the inclusion of the 
deformed shape of the nucleus in the calculations. One can see that the 
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Nucleus 

Po 

Pi 

P2 

P3 

P4 

P5 

P6 

P7 

Ps 

P9 

spherical 

0.04 

0.03 

0.01 

0.01 

0.43 

0.31 

0.62 

0.31 

0.71 

1.00 

deformed 

0.14 

0.48 

0.62 

0.32 

0.56 

0.14 

0.22 

0.18 

0.16 

1.00 


Nucleus 

Pio 

Pii 

P 12 

Pl3 

Pl4 

Pl5 

P 16 

Pl7 

P 18 

spherical 

1.00 

1.00 

1.00 

0.67 

0.09 

0.01 

0.91 

1.00 

0.01 

deformed 

1.00 

0.72 

0.38 

0.02 

0.23 

0.51 

1.00 

1.00 

0.01 


Table F.9: The fusion probabilities for the capture of the a-particle by the spherical 
and deformed nuclei “^^Ca obtained by the MIR approach (parameters of calculations: 
10000 intervals at rmax = 70 fm; errors in obtaining results for the deformed nucleus are 
ei = 0.0191, £2 = 0.0263 and £3 = 0.0040). 


MIR approach allows us to obtain such coefficients for the deformed nuclei 
with an accuracy similar to that of the spherical nuclei. The corresponding 
cross-section for the a-capture by the deformed nucleus ^^Ca with included 
fusion probabilities is given in Fig. lF.13l One can see that line 3, obtained for 
the deformed nucleus describes the experimental data with good accuracy, 
similar to that of line 2 obtained for the spherical nucleus. We observe the 
presence of small oscillations in the calculated spectrum after the inclusion 
of the nuclear deformation. 
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Figure F.13: The capture cross-sections of the a-particle by the ^“^Ca nucleus 
obtained by the MIR method (parameters of calculations: 10000 intervals at 
r-max = 70 fm, parametrization H). Here, the data labeled 1 are the experimental 
data extracted from [H^, solid blue line 2 is the cross-section at / ma x = 17 for the 
spherical nucleus, dash-dotted purple line 3 is the cross-section at /max = 17 for the 
deformed nucleus, and dash dotted green line 4 is the cross-section at /max = 15 
for the spherical nucleus (cross-sections are defined by (IF.lh . where /max is the up¬ 
per limit in the summation). Lines 2 and 3 are obtained with the included fusion 
probabilities, line 4 without the fusion probabilities. 
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